Non-contact gears: I. Next-to-leading order contribution to lateral Casimir force 

between corrugated parallel plates. 
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We calculate the lateral Casimir force between corrugated parallel plates, described by ^-function 
potentials, interacting through a scalar field, using the multiple scattering formalism. The contribu- 
tions to the Casimir energy due to uncorrugated parallel plates is treated as a background from the 
outset. We derive the leading- and next-to-leading-order contribution to the lateral Casimir force 
for the case when the corrugation amplitudes are small in comparison to corrugation wavelengths. 
We present explicit results in terms of finite integrals for the case of the Dirichlet limit, and exact 
results for the weak-coupling limit, for the leading- and next-to-leading-orders. The correction due 
to the next-to-leading contribution is significant. In the weak coupling limit we calculate the lateral 
Casimir force exactly in terms of a single integral which we evaluate numerically. Exact results for 
the case of the weak limit allows us to estimate the error in the perturbative results. We show that 
the error in the lateral Casimir force, in the weak coupling limit, when the next-to-leading order 
■ contribution is included is remarkably low when the corrugation amplitudes are small in compari- 

son to corrugation wavelengths. We expect similar conclusions to hold for the Dirichlet case. The 
analogous calculation for the electromagnetic case should reduce the theoretical error sufficiently for 
Q ■ comparison with the experiments. 

CD 
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(N ' I. INTRODUCTION 

> 

\ The Casimir force, as exhibited between neutral metallic parallel plates, was discovered theoretically in 1948 [1]. 
The Casimir torque between asymmetric materials was first studied in 1973 [2]. Recently, theoretical study of the 
lateral Casimir force between corrugated parallel plates was pioneered and developed by the MIT group in [3-7]. In 

ly-j , particular, in [6], the authors evaluated analytic expressions for the lateral Casimir force, to the leading order, between 
two corrugated parallel plates perturbatively. Experimentally, the Casimir interaction between corrugated surfaces 
was explored during the same period of time by Roy and Mohideen in [8]. This experiment measured the lateral 
Casimir force between a plate, with small sinusoidal corrugations, and a large sphere with identical corrugations. The 
motivation was to study the nontrivial boundary dependence in the Casimir force. The experimental parameters in our 

■ ■ »h , notation are (see figure 1): h — 60 nm, d = 1.1 /im, and a = 0.1 — 0.9 /im, where h is the height of the corrugations, d 
is the wavelength of the corrugations, and a is the mean distance between the plates. The corresponding dimensionless 
quantities are: k^a = 0.6 — 5.1, ^ = 0.07 — 0.6, and kgh = 0.3, where fco is the wavenumber related to the spatial 
- ' wavelength of the corrugations. 

Experimental data was analyzed based on the theoretical results obtained from the proximity force approximation 
(PFA), and has been presented in [9, 10]. The validity of the PFA in the above analysis has been the topic of a 
recent debate and controversy, see [11—13]. Theoretical results based on perturbative approximations as done in [5] 
do not settle the issue because the error keeping only the leading order may be high. It is generally believed that the 
next-to-leading-order calculation will be able to throw light on the issue. We carry out this calculation for the case 
of scalar fields. The analogous calculation for the electromagnetic case should now be straightforward. 

This paper in principle is an offshoot of [14] where we shall deal with corrugated cylinders to study non-contact 
gears. While evaluating the leading order for the case of corrugated cylinders it was noticed that it would be possible to 
extend the calculation to the next-to-leading order. This led to the study in the present paper. In this installment we 
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present the next-to- leading-order calculation for the case of corrugated parallel plates. The leading order calculation 
for the corrugated cylinders, which in itself is a significant result, will form the sequel [14] of this paper. The next-to- 
leading-order calculation for the corrugated cylinders is in progress. 



II. FORMALISM 



In this section we shall describe the formalism and derive the key formula used for calculating the Casimir energy. 
This has been done in various papers before, (see [15], [16], [17], and references therein). We hope our derivation 
using Schwinger's quantum action principle techniques will be illuminating. In an earlier paper [17] describing the 
multiple scattering formalism it was mentioned that the use of the scattering matrix, T, was equivalent to using the 
full Green's function, and required the same computational effort. As a justification of this comment we exclusively 
use the full Green's function in this article. 



A. Vacuum energy in field theory 



Let us consider a scalar field, (j)(x), interacting with a scalar background potential, V(x), described by the Lagrangian 
density 



=~\ d^{x)d^{x) - \ V(x) <t>{x)\ 



(1) 



In terms of the source function, K(x), corresponding to the scalar field, we write the action for this description to be 

W[(f>;K} = j d 4 x K(x)<j)(x) +C((j>{x)) . (2) 
The vacuum to vacuum persistence amplitude, 

Z[K] = (0+|0_>* (3) 
which generates the amplitudes for all the physical processes, satisfies Schwinger's quantum action principle, 

SZ[K] = i({)+\5W[<t>:K] \Q-) K . (4) 
Our immediate task will be to get a formal solution for the vacuum amplitude, Z[K], in the form 



Z[K] 



JW[K] 



(5) 



where W[-K], which is not an action anymore, is dependent only on the source function. Note that the action, W[(j>; K] 
in cq. (2), which satisfied the action principle, was described in terms of both the (operator) scalar field and the source 
function. 

Variation with respect to the source function in the quantum action principle in eq. (4) allows us to write 



= (0 + |^)|0-)* 
^ [ > ~ Z[K] 



1 1 6Z[K] 



Z[K] i 6K(x) 



(6) 



where the redefined scalar field, on the left of the above expression, is an effective field. This can be used to replace 
operator fields with functional derivatives with respect to the sources. Variation with respect to the scalar field in 
cq. (4) gives us 



d 2 - V{x) 



1 6Z[K] 



K(x)Z[K] 



i SK(x) 

which can be immediately inverted and written in the form, after using eq. (6), 

1 1 5Z[K] 



ip(x) = 



Z[K] i 6K(. 



Q= [ d 4 x'G(x,x')K(x'), 
x) J 



(7) 



(8) 
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where we defined the inverse of the differential operator, the Green's function, as 

- 'd 2 -V(x)] G{x,x')=5^\x-x'). 



(9) 



The solution to eq. (8) is a Gaussian in the functional sense, and when written in the desired form in cq. (5), lets us 
identify 



W[K] = Q[V] + \j dAx j d 4 x'K(x)G(x,x')K(x'), 



(10) 



where Q[V] is a functional of the background potential alone. For the case when the background potential is switched 
off, described by V = 0, we have Z [K] = cxp(iW [K]), where 



W [K] = Q[0] + d 4 x J d 4 x'K(x)G {x,x')K(x') 



(11) 



where, Gq{x,x') is the corresponding Green's function in eq. (9) for the case when V = 0. Now, in the absence of a 
source function the vacuum should not decay, which amounts to the statement (0 + |0_) x=0 = 1. This implies that 
Wo[0] = which when used in cq. (11) lets us conclude that Q[0] — 0. 
Variation with respect to the background potential in cq. (4) yields 

5Z[K] = - l -J d 4 x SV(x) (0+ 1 <f>(x)<f>(x) |0_ ), (12) 

where we can replace the operator fields with functional derivatives and then write the formal solution as 

Z[K] = e -5 /d4xy(x) 7^TTOW e 4/^/rfVK(x)Go(x ) x')X(x') i ( 13 ) 

Using the standard identity, we have 

W\K] = In GGq 1 + \j dAx j d 4 x'K(x)G(x,x')K(x'), (14) 

which uses (1 + GoV)~ 1 Go = G. Using eq. (14) in eq. (5) we observe that in the presence of a background the vacuum 
to vacuum transition amplitude is not unity, instead it evaluates to 



<0+|0_> 



K=0 



iW[0] 



e |Tr In GGo\ 



(15) 



For the case when the process under investigation is time independent, denoting r to be the time associated with 
the physical process, we can evolve the vacuum state using the Hamiltonian of the system and thus conclude 



(0+|0_) K=u = (0_|e-^ T |0_) =e 



-iEt 



(16) 



where we assumed the vacuum to be an eigenstate of the Hamiltonian. Comparing the two forms, eq. (15) and eq. (16), 
we thus identify the energy of the vacuum in the presence of a background to be 



E= — Tr In GG^ 1 
2r 



(17) 



which serves as the central formula for calculating the Casimir energy. Further, for the case of a time independent 
situation, making a Fourier transformation in the time variable and using translational independence in time, we can 
formally write 



E 



i r du 
~2J 2^ 



Trln 



- lu 2 - V 2 + V 



+ 



1 f du) 

2 J 2^ 



Trln 



- u? - V 2 



where we used eq. (9) to define the Fourier transformed Green's function in the form G 1 = [- 
Integrating by parts and throwing away the boundary terms we derive 



E 



— 2lj 2 Tr G + - 
2tt 2 



dui , 
2^ 



2w 2 TrG , 



(18) 

-n 



(19) 



which is an alternate expression for evaluating the Casimir energy. This latter form has been used in studies [18] 
related to surface divergences which requires the evaluation of the energy density rather than the total energy. 
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B. Statement of the problem 



We consider two semitransparent, corrugated plates, parallel to the x — y plane, described by the potentials, 

V i (z,y) = X i 8(z-a i -h i {y)), (20) 

where i = 1,2, are labels that identify the individual plates, and we define a = a 2 — ai > 0. The functions hi(y) are 
designed to describe the corrugations associated with the individual plates. Let us define the function 

a(y) = a + h 2 (y)-h 1 (y) (21) 

which measures the relative corrugations between the two plates. We shall define the corrugations hi{y) so that the 
mean of the relative corrugations is a. In general we require a(y) > 0, which is a restriction on the corrugations. 
Translational invariance is assumed in the x direction. 

The change in energy due to the change in the mean distance between the plates, a, leads to the conventional 
Casimir force which points in the direction perpendicular to the surface of the plates, and is expressed as 

F Cas = --, (22) 

where E is the total Casimir energy associated with the corrugated plates including the divergent contributions 
associated with the single plates. The divergent contributions being independent of a, do not contribute to the 
Casimir force. When the plates have corrugations on them we expect to have a change in the total energy due to 
a shift of one of the plates parallel (lateral) to the other plate. The force corresponding to the change in the total 
energy due to this shift acts in the lateral direction and is called the lateral Casimir force. The shift is mathematically 
described by a translation in the y-coordinate, h\{y + yo), which corresponds to a lateral shift of one plate with respect 
to the other. The lateral force is expressed as 

FlE 

^Lat = -%-, (23) 
dyo 

where yo is the measure of the translational shift. 

We note that there will be no lateral force between the plates if the corrugations are switched off by setting hi(y) = 0, 
i = l and 2. The physical quantities associated with this configuration will thus act as a background, and a reference, 
and we shall find it convenient to denote them by the superscript (0) to mean zeroth order. The potential for the 
background will thus be described by 

V^\z) = \ l d(z-a l ), (24) 

which has no y dependence. The total Casimir energy associated with the background, due to the two uncorrugated 
plates, will be denoted as E^ , which will include the divergent contributions associated with the single plates. This 
background energy will be independent of the displacement yo due to the y independence of this configuration. Thus, 
we can conclude that 

We further note that there will be no lateral force between the plates if either one of the plates have their corrugations 
switched off by setting hi(y) = 0, i = 1 or 2. The Casimir energy associated with this configuration can be written as 

Ei = E- E (0) (26) 

where Ei is the additional contribution to the energy with respect to the background energy due to the presence of the 
corrugations on one of the plates. Throughout this article we shall use A to mean the deviation from the background. 
For example, we will have AVi(z,y) = Vi(z,y) — V^°\z). Using the argument of y independence, or translational 
symmetry, we can again conclude that 

BE- 

(F Lat ) l = --^=0. (27) 
oyo 

In light of the above observations we are led to write the total Casimir energy, for the case when both the plates 
have their corrugations switched on, in the form 

AE = E- E (a \a) = E^a, h{) + E 2 {a, h 2 ) + E 12 (a, h^h*, yo), (28) 
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where E\2 is the contribution to the total energy due to the interaction between the corrugations in the plates. Only 
this part of the total energy contributes to the lateral Casimir force. Thus, we conclude that 

F Lat = -^A£=-^. (29) 

Our central problem will be to evaluate E12 using eq. (17). 



C. Casimir energy contributing to the lateral force 

Using the central formula derived in the multiple scattering formalism in cq. (17) to evaluate eq. (28) we have 
AE = E- E<® = — TrlnGGrT 1 - — TrlnG^Gn 1 = — Tr In GG^~\ 



2r 



2r 



2t 



(30) 



where Go is the free Green's function introduced in eq. (11). Note that Go cancels in the above expression and the 
reference is now with respect to the uncorrugated plates. The Green's function G satisfies eq. (9) with potentials 
described by eq. (20), which in symbolic notation reads 



-d 2 + V 1 +V 2 



G = l, 



and the corresponding Green's function associated with the background satisfies the differential equation, 

-a 2 + y 1 (0) + T/ 2 (0) ]G(°) = i. 

The above two equations can be used to deduce 



GG (o)- 



l + G (0) AUi + G (0) AU 2 



Next, one makes the observation that the above expression can be rewritten in the form 

GG (orl = G 2 G (orl [l - G 1 AV 1 G 2 AV 2 \ ^GiG^" 1 , 



(31) 



(32) 



(33) 



(34) 



where Gj (i = 1, 2) are the Green's functions for the parallel plates when only one of the plates has corrugations on 
it. The differential equations for G»'s are 



-d 2 + + U 2 (0) + AvJ Gi = l, 



(35) 



which together with eq. (32) can be used to deduce 

G,G (orl = (1 + G^AVi)- 1 = 1 - G (0) AU,(1 + G^AV;)" 1 = 1 - (1 + G (0) AV i )- 1 G {0) AV t . (36) 

Using eq. (34) in eq. (30) we we immediately obtain eq. (28). The last term in eq. (28), E12, which is the only term 
that contributes to the lateral force between the two corrugated plates, can be read out from eq. (34) to be given by 
the expression 



E 12 = -T^Trln 
It 



1-G 1 AV 1 G 2 AV 2 



(37) 



We could have written this down immediately, but we hope it was instructive to go through the steps because it clarifies 
our notation, which would anyhow have required us to write many of the above equations, and also to emphasize that 
the configuration due to parallel plates is treated as a background from the outset. 
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D. Formal series expansion 



Formally expanding the logarithm in the above expression and using eq. (36) to expand GjAVj in terms of G^AVi 
we can write 



i 1 oc oo 

^ = 27^E^ E E (-iTU-!)™ 2 {g (0) aVi j™ 1 |g (0) AV^j' 



m — 1 ]_ m =1 712=1 

Our potentials in eq. (20) can be formally expanded in powers of hi in the form 



A^) = 5>^(z, y ) = E 



-My)]" d n T/(0) 



n— 1 n— 1 

so we can further write the series expansion as 



n! <9z 



e -fc 4 (»)& _ i 



<5(z - a*), 



i - 1 OO OG 

B « = ^Tr £ 1 £ £ (_ Al )" 1 (-A 2 )- {g<°> [e^ 8 - 1]^}" {g<°> [e"^ 9 - l]<5 2 }" 

m—1 L^i— 1^2 — 1 

where we again use symbolic notation and suppress the variable dependence in hi, <9, and the delta functions. 



(38) 



(39) 



(40) 



III. LEADING ORDER CONTRIBUTION 

In this section we shall derive the leading order contribution to the Casimir energy when the corrugation amplitude 
is small in comparison to the corrugation wavelength. This has been evaluated for the Dirichlet scalar case and 
the electromagnetic case in [5, 6]. We obtain the result for scalar S- function potentials as a warm up exercise in 
preparation for the calculation for the next-to-leading order. In this section we shall illustrate our methodology which 
will be further used in the higher order calculation. 



A. Second order perturbation in Casimir energy due to corrugations 

For the particular case when the corrugations can be treated as small perturbations we can approximate the 
potentials by keeping a few terms in the expansion in eq. (38) where we use the superscripts (n) to represent the n-th 
order perturbation in a quantity. Thus, to the leading order the interaction energy of the corrugations in eq. (38) 
takes the form 



E® = -Tr 
12 2r 



G<°>AVi (1) G< >AV 2 (1 



(41) 



We observe that the potentials in the n-th order are given by derivatives acting on <5-functions in eq. (39) which can 
be transferred over to the Green's functions after integration by parts. The background Green's function, which is a 
solution to eq. (32), can be written as 



G(°)(x,i,x',i')= / — e — (*-*') / ^i e *M*-*') / ^ e *M»-*V 0) (*,«';«), 
J 2-k J 2tt J 2tt 



(42) 



where n 2 — k 2 + k 2 , with k 2 = k 2 — uj 2 , where k x . y and lo are the Fourier variables corresponding to the space-time 
coordinates x, y and t respectively. The reduced Green's function, g^°\z, z'\ n), satisfies the equation 



- k 2 - Xi5(z - oi) - X 2 5(z - o 2 ) 



g^(z lZ '-K) = 8(z-z'). 



In terms of the reduced Green's function defined above we can write eq. (41) as 



E 



(2) 
12 



= r r°° ^~ hi{ki _ k2) ~ h2{k2 _ ki)L {2 ){kM 



(43) 



(44) 
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where L x is a large length in the x direction and where hi{k) are the Fourier transforms of the functions hi(y), which 
describe the corrugations on the parallel plates, 



/oo 
dye~ iky hi(y). 
-oo 

The kernel L^ 2 \ki, k 2 ) in cq. (44) is suitably expressed in the form 

L { - 2 \k u k 2 ) = - — iidii I {2 \ Ki ,k 2 ), 
47r Jo 

where we have switched to imaginary frequencies by a Euclidean rotation, w — > i(, and defined 

d d 



I^( Kl , k 2 ) = AxA 2 |- ^--[9 {0) {z, z; Kl ) ff (0) (z, z; k 2 ) 



(45) 



(46) 



(47) 



z— ai ,z— a,2 



where «? = k 2 + kf. Using the reciprocal symmetry in the Green's function in cq. (47) we deduce that 

I {2 \k 1 ,k 2 ) = I { - 2 \k 2 ,k 1 ). 

We evaluate the derivatives in the expression for I^(m, k 2 ), using the prescription described in appendix A, as 



(48) 



/ (2) (*i,k 2 ) 



Ai A 2 e - a < Kl+K2 ) 
2m 2n 2 AiA 2 



Ai \ / A 2 \ o Ai \ /. A 2 



where Aj's are given by eq. (A2) after replacing n — > m. 



(49) 



Dirichlet limit 



For the case of the Dirichlet limit (aAi, 2 3> 1) the expression for I^(ki,k 2 ) in eq. (49) takes on the relatively 
simple form 



I^\k u k 2 ) = -- 



K 2 



sinh ma sinh n 2 a 

where the subscript D stands for Dirichlet limit. Using the above expression in eq. (46) we have 



L ( £ > (k 1 ,k2) = [ ndn^p- r 

4-7T J sinh ki a si 



K2 

sinh K 2 a 



(50) 



(51) 



We describe the correspondence of the expression for the interaction energy, in eq. (44), in the Dirichlet limit, to the 
result in [6] in appendix B. 

We shall digress a little to answer how the above result is related to the Dirichlet Green's function. We start by 
recalling the reduced Green's function corresponding to eq. (43) in the Dirichlet limit, which can be derived by taking 
the Ai i2 — > oo limit in g 2 ) \z, z'; k) in cq. (Alb) defined in region 2 (ai < z, z' < a 2 ) in figure 10, 



(°)/ ' \ 
9 D (z,z ;k) 



sinhK(z< — ai) sinhK(a 2 — z>) 



(52) 



K sinh na 

where z < and z > stand for Min(z, z') and Max(z, z') respectively. It is trivially verified that the above function 
satisfies the Dirichlet boundary conditions, g^,\a{, z'\ k) = and g^\z, k) = 0. Using eq. (52) we can evaluate 

d d 



dz<- dz-~ 



sinh K,a 



(53) 



J z < =a 1 ,z > =a 2 

Using the above result eq. (50) can be expressed as the product of derivatives of two Dirichlet Green's functions. It 
is also worth mentioning that the result in cq. (50) can also be derived by exclusively using g^{z,z']K) in eq. (47) 
after taking the limit Ai, 2 — * oo. This is expected because the result in the Dirichlet limit should depend only on the 
quantities between the plates. Thus, all the results related to the Dirichlet limit can be derived without relying on 
the averaging prescription for taking derivatives described in appendix A. However, for the more general case being 
considered here, we require the averaging prescription to derive eq. (49). 
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FIG. 1: Parallel plates with sinusoidal corrugations. 



2. Weak coupling limit 

For the weak coupling limit (aAi.2 <C 1) the expression for I^ 2 \k\, K2) in eq. (49) takes the very simple form 



(54) 



where the subscript W stands for the weak limit. We point out that in the case of weak coupling the averaging 
prescription for the evaluation of the /-kernels was not necessary because only the free Green's functions come in. 
Using the above expression in eq. (46) we can evaluate the corresponding L-kernel as 



Ly(ki,k 2 ) = — — 

16tt Jq 



RdR ( Kl+K2 ) 2 e -a( Kl + K2 ) = l_ 

K1K2 167T da 



i e -o(|fci|+|fe 2 |) 



(55) 



where in the evaluation of the integral we used the change of variables n± + k 2 = x, and the corresponding relation 
R, dR (ki + k 2 ) = K1K2 dx. 



B. Sinusoidal corrugations: Leading order 

For the particular case of sinusoidal corrugations, as described in figure 1, we will have 

hi(y) = hi8in[ko{y + vo)], 



(56a) 
(56b) 



where ko — 2-k /d is the wavenumber corresponding to the corrugation wavelength d. We get nonzero contributions, to 
the leading order, only for the case when both plates have the same corrugation wavelength. The Fourier transforms 
hi(k) for sinusoidal corrugations get contributions in the form of delta functions 



h^k) = hx^r \e ikoyo 5(k - k ) - e- lkoVo S(k + ho) 
2i L 



(57) 



Using the above expression in eq. (44), and after interpreting 2ttS(0) — L y to be the infinite length in the y direction, 
we write 



, h\h 2 
= cos k y — — 

47T 



L x L y 



dkL { - 2) (k,k + ) = -cos/cq 



Vo 



hih 2 
16^2 



dk 



k dn 1^ (k, k+), 



(58) 



where we have used the symmetry property in the /( 2 )-kernel, noted in eq. (48), and performed suitable rescaling in 
the integration variables. We have used the notations k± — k ± fc , and n 2 ± — R 2 + k±. The I^ 2 ' -kernel in the above 
expression is provided by eq. (49). 
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FIG. 2: Dirichlet limit: Plot ol A^'Qtoa) versus koa 



1. Dirichlet limit 

In the Dirichlet limit, where the /( 2 '-kernel is given by eq. (50), the interaction energy in eq. (58) can be expressed 
in the form 



E i2 n . „. \ 71-2 h i h i 



L X Ly 



cos(fc j/o) 



240 a 3 a a D 



(59) 



where the function A^ ) ,1 \koa) is normalized such that A^'^O) = 1, which becomes a convenience when we compare 
the results with those obtained by the proximity force approximation. The expression for A^'^(fc a) is 



4^(fc a) = ^|- dkL%\k,k 

^ J — OO 



1 60a 5 ,-„/-_ ,__ (2) , 



4V 7T" 



— oo JO 



dk / ndnI D (k,k + ) 



(60) 



(2) (2) 

where the -kernel, and the corresponding L D -kernel, in the Dirichlet limit were given in eq. (50) and eq. (51) 
respectively. Using the notations, s 2 — s 2 + 1 2 , s\—s 2j r{t± t ) 2 , and t = k a, we can write 



a™ to) 



15 



.10 



s ds dt 



sinh s sinh s + 



(61) 



This expression can be numerically evaluated and has been plotted in figure 2. The value of the function at to = 
can be evaluated exactly by rewriting the integral in terms of spherical polar coordinates to yield 



4^(o) 



30 



poo 

/ s 2 ds 


s 


10 


sinh s 



(62) 



It should be mentioned that the result corresponding to eq. (61) for the electromagnetic case and the scalar Dirichlet 
case have been evaluated exactly in [6]. 

The lateral Casimir force can be evaluated using eq. (59) in the definition of lateral force in eq. (29) which yields 



( 2) 

^Lat.D = 2k o a sm ( fc o2/o) 



F 



(0) 



Cas,D 



hi h 2 .(1,1)/, v 
~~a~a D 



(63) 



where l-F^g D | is the magnitude of the normal scalar Casimir force between two uncorrugated parallel Dirichlet plates 
given as 



p(o) 

r Cas,D 

L X Ly 



480 a 4 ' 



(64) 
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FIG. 3: Weak coupling limit: Plot of A<$ (k a) = e- k ° a e 2 {k a) versus koa. 



2. Weak coupling limit 



Let us start by recalling the expression for the Casimir energy and Casimir force between uncorrugated parallel 
plates in the weak limit: 



E, 



(o) 

Cas,W 



L X Ly 



A X A 2 1 
~32tt 2 a 



and 



r Cas,W 

L X Ly 



AiA 2 1 
32 7T 2 a 2 



(65) 



In the weak limit we have seen that it is possible to evaluate the L- kernel in eq. (55) without any effort. Using this 
in eq. (58) we have 



E 



(2) 

12,1V 



L X Ly 



cos(k Q y 



AiA 2 hih 2 d 
16 7r 4 7r da 



dk e ~ a (\ k \+\ k + k °\) 



which after performing the integral can be written in the form 



( 2) 

^12, w = cos ( k oyo) 

where we have introduced the function 



E, 



(0) 

Cas.W 



h\ h 2 (2) , . 
a a 



A2) { ,4 d 2 



1 

— e 
to 



J.IIL 



e-2{to) 



(66) 



(67) 



(68) 



where e n {x) = X)m=o ^TT ^ s ^ nc ^ ne truncated exponential function which approximates to unity for x <C n. A^) (to) 
has been plotted in figure 3. Using eq. (29) we evaluate the lateral Casimir force in this perturbation order as 



(2) 

^Lat.W = k o a sin(/c ?/o) 



F, 



(0) 



Cas.W 



hi h 2 ( 2 ) n s 

A W ( k o a ) 

a a 



(69) 



IV. NEXT-TO-LEADING ORDER CONTRIBUTION 



We start by collecting the terms in eq. (38) in the form 



El2 — &12 



E 



(3) 



+0{h/af. 



J 12 ' ^12 

To economize typographical space and for bookkeeping purposes we introduce the notation 



[G<°>AV]g 



Tr 



G (0) AU (mi) G (0) AU (m2) • • • G (0) AU (mjv) 



(70) 



(71) 



11 



where the superscripts m n denote the m n -th order contribution in the power series expansion of the potential as 
defined in eq. (39). The subscripts i n identifies the potential that is contributing. In the case under consideration we 
have only two potentials coupling and thus i n = 1, 2 for any n. As an illustrative example we have 

[G (0) AVlg;ig = tAg^av^g^av^g^av^] . (72) 

Using the above notation we can write 

E[f = A{ [G (0) Ay] (V) + [G (0)AV]jl : »j - [G< 0) AV]jl;$ - [G^AV]^}. (73) 

We argue that for the case of sinusoidal corrugations of the same wavelength on both plates the third order 
perturbation does not contribute. (Contributions from third order are nonzero if the wavelength of the corrugations 
of one plate is double that of the other plate.) Since the plates under consideration arc infinite in extent, the interaction 
energy is independent of translations in the y direction. Let us now make a translation of the amount koy = ir in the 
potentials in cq. (56). This amounts to replacing hi — > —hi. Thus invariance under translation requires that the total 
power of hi 8 in the result should be even. This rules out the third order perturbation from contributing. 



A. Fourth order perturbation in Casimir energy due to corrugations 

Starting from eq. (38) and keeping terms contributing to the fourth order we have 

E^=E^+E[^+E[r\ (74) 

where the superscripts represent the powers of /ij's. For example, the superscript (3, 1) represents terms involving 
h\h 2 - All terms except one gets contribution from m = 1 in the logarithm expansion in eq. (38). Since the m = 2 
term is distinct in structure we further separate it out by defining 

e[^=e[^ A + e[^ B , (75) 

where the term involving the superscript A is the contribution from m — 2 in (38). Using the notation in eq. (71) we 
can thus collect the terms contributing to eq. (74) and eq. (75) as 



ci(2,2)A 
^12 




}■ 






(76a) 


^12 


= A{ [G (0 Ay] (M ; M)_ 


[G<°>AV]g;$ 


- [G^AV]g : ig 


+ [G (0) AV]g3} > 


(76b) 


£(3,1) 
^12 




[G<°>AV]g : $ 


- [G (0) AV]g;$ 


+ [G (0) AV]g;iJ} > 


(76c) 


^12 


= ^MSS- 


[G<°>AV]g : $ 


- [G^AV]g : ig 


+ [G< )AV]g;g}. 


(76d) 



In terms of the reduced Green's function defined in eq. (42) and eq. (43), and the Fourier transform of the corrugations 
in eq. (45), we can write each of the terms in eq. (76) in terms of a corresponding L- kernel. In terms of the notation 
introduced in eq. (71) these will read as 

( _i)" jl [G (o) Av] = Lx j^...^~ h T:(h-k 2 )... k: (k N - h) L^-^Hh, . . . , kN ), (77) 

where implicitly we have interpreted the powers of the Fourier transformed corrugations as 

h^(h-k 2 )= J ^...^^h ll (k 1 ~ qi )h ll (q 1 -q 2 )---h n (q mi ^-k 2 ). (78) 

Observe that each term contributing to the M-th order (M = mi +m 2 + . . . + mjv) has exactly M Zi's. Also, the total 
number of variables inside the L-kernel is N which is in general less than or equal to M. It might be appropriate to 
call them Appoint kernels. Proceeding further, in the spirit of the second order calculation in eq. (46), we introduce 
the corresponding /-kernels for the L-kernels as 

L fcZM N) ^> ■- kN) =--kl RdR fci*7V. • ■ • - (79) 
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where the /-kernels are expressed in terms of derivatives operating on the reduced Green's function defined in eq. (43) 
as 



j(mi,...,m») 
(il ,...,ijv) 



(«i, ...,k n ) = (-1) 



A (9 mi 



N ^ii _ 

mi! mAr! dz™ 



Qm N 

uz N 



g {a) (z N ,zi,Ki) ■ ■■g w (z N -i,z N ;K N ) 



(0), 



(80) 



For clarification wc illustrate the evaluation of a particular term which should also serve as an illustration of the 
notation. The term we consider is 



A [G (0) AV1 (W) 



J (l,l,2) 



L, 



'^k'^k'^K ~ ^^lfa ~ fc 3)ft-2(fc3 - illijlj^l, ^2,^3) 

dh dh dh dh _ _ ^/^^ _ ^-j ^ 2 ( fc4 _ fcj) L^' 2 '^(fc 1; fc 2; fc 4 ). (81) 

Z7T Z7T Z7T Z7T V , , / 



Notice how the particular k dependence in the L-kernel is unambiguously specified. The corresponding /-kernel using 
eq. (79) and eq. (80) is given as 



7 (i'i,2) K 2, K4) = --^^7^2 ^7 f5 (0) (^3, zi ; ki) 3 (0) (zi, z 2 ; K2)g (0) ( z 2, z 3 ; k 4 ) 



2 <9zi <9zf 



(82) 



zi=a 1 ,z 2 =a 1 ,z 3 =a 2 



Using the reciprocal symmetry in the Green's functions we can learn the following symmetries in the /-kernels 
associated with the terms in eq. (76): 



r(i, 1,1,1)/ \ r (i, 1,1,1)/ \ 

^(1,2,1,2) ( Kl ' K2 ' K3 ' K4 ) = ^(1,2,1,2) V K 2, Kl, «4, 

r(l, 1,1,1)/ \ _ r(l, 1,1,1)/ \ 

-'(1,1.2,2) V Kl ' K 2' K 3, K4) — ^/ lil;2i 2) l K 3, «2, K\, K4), 

r(l,l,l,l)/„ ,„ „ ,„ x, _ r (l,l,l,l) 
(1,2,2,2) 

r(l,l,l,l)/ 



r(l,l,2)/ \ r(l,l,2)/ \ 

J (l,l,2)( Kl ' K2 ' K 3) = J( 1)li2 )(K3,K2,Kl), 



r(2,l,l)/ \ \ 

(1,2,2) l K l, K 2, «3j — -'(1,2,2) \ K 2,Kl,K 3 ) 



(2,1,1), 



T^±,-L, ±,-L_) / \ T-l^i, 1,±, 1 J / \ 

(1 2,2 2) V^l' ^4 J — 2 2 2) v^ 2 ' 1' 4' 



/, 



(1.2.1) / \ \ 

(1.2.2) l K l, K 2, «3j — -'(1,2.2) V K 2, «1, «3j, 

r(l,2,l)/ \ r(2,l,l)/ \ 

J (1,1,2)( K 1' K 2,K3) = J(i,i 1 2)( K 3,«2,Kl)- 



_ r(l,l,2). 



tII.1,1.11/ \ x(l, 1,1,1)/ \ 

/ (l,l,l,2)" Kl ' K2 ' K 3: «4) = ^(1,1,1,2) ( K 4, «3, «2, 



(83a) 
(83b) 
(83c) 
(83d) 



We have not listed the 2-point /-kernels in the above list because they have the symmetry of the kind in eq. (48). In 
fact any 2-point kernel will have the following symmetry: 



7 -(mi,m 2 )/ \ _ j\mi,m2)( \ 

■'(1,2) \ K ^ K 2) — -'(1,2) l K 2,«lJ- 



(mi,m 2 ) , 



(84) 



In the above discussion involving very general notations we have in principle expressed each term in eq. (76). 
Specific evaluation of each term involves the evaluation of the corresponding /-kernels which are given in terms of the 
derivatives of the Green's functions. The derivatives are evaluated using the prescription described in appendix A. 
We can thus collect the terms in eq. (76) into four distinct L-kernels in the form 



E 



(2,2)A 
12 



E 



(2,2) B 
12 



dhdk^dk^dk^^ _ k 2 )h 2 (k 2 - k 3 )hi(h - h)h 2 (h ~ ki) L (2 ' 2)A (fci, k 2 , k 3 , fc 4 ), 

Z7T Z7T Z7T Z7T 



E 



(3,1) 
12 



r dh 
J ~^ 

r dh 

J 27 



ak^dk^dh _ fej^fe _ k 3 )h 2 (k 3 - ki)h 2 (ki - h) L {2 - 2)B (h, k 2 , k 3 , fe 4 ), 

Z7T Z7T Z7T 

ak^dk^dh _ _ ^^(fcg _ h)h 2 (h - h) i (3,1) (fci, fc 2 , fe 3 ,fc 4 ), 

Z7T Z7T Z7T 



(1,3) 
12 



rffci dfc 2 dfc3 dfc4 ~ 



2tt 2tt 2tt 2tt 



fti(fei - k 2 )h 2 (k 2 - k 3 )h 2 (k 3 - k 4 )h 2 (k 4 - h) L {1 ^\h, k 2 , k 3 , fc 4 ), 



(85a) 
(85b) 
(85c) 
(85d) 



where it should be noted that different L's combine with specific combination of h's. The factor of one-half in eq. (85a) 
can be traced back to the coefficient of the second term in the expansion of logarithm in eq. (38). The respective 
kernels L^" 1 ^ above are related to their corresponding /( m ' n ) by eq. (79), which are given by 

/ (2 ' 2)a (ki, k 2 , k 3 , k 4 ) = I^I'1'^(ki,k 2 ,K3,K4), (86a) 

/ (2 ' 2)B (ftl, K 2 , K 3 , Ki) = ^(1,1'^) (Kl, «2, «3j «4) + 4l'l,2) ( K l> K 2, K3) + ^(1,2,2) ( Kl ' «3, «4) + ^(1,2) ( K li ( 86b ) 

r(3,l)/ \ r(l, 1,1,1)/ \ i r(l,2,l)/ \ . r (2,l,l)/ x , r (3,l)/ n , Qa x 

I y ' '(Kl, K 2 , K 3 , Ki) = I( 1112 ){Ki,K 2 ,K 3 ,Ki) + i ( \ )1|2 ) («1, «2, K4) + i ( \ |1)2 ) («1, K3, «4) + i( 1)2 ) («1, K4), (86c) 

r(l,3)/ \ r(l, 1,1,1)/ \ 1 r(l,2,l)/ \ . r (l,l,2)/ n . r (l,3)/ x 

/ l ' '(Ki, K 2 , «3, K4) = ^(1,2,2,2) ( K 1' K 2,K3,«4) + ■f( 1)2l 2) K 2, K4) + i ( \ |2)2 ) («1, «2, K3) + ■'(1,2) ( K l> K 2)- (86d) 
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Proceeding in the spirit of section III has brought us to the point of evaluation of the thirteen /-kernels on the 
right hand side of eq. (86). Using the prescription described in appendix A we have evaluated all the thirteen 
kernels and explicit expressions analogous to eq. (49) have been derived. This would have involved a lot more labor 
and bookkeeping if not for the facilitation achieved by the use of Mathematica. We shall not display the explicit 
expressions here because they are too long. But, as in section III, these expressions simplify considerably in the 
Dirichlet limit and the weak coupling limit. 



1. Dirichlet limit 



Observe that the reduced Green's function defined by eq. (43) has a well defined Dirichlet limit. In the light of this 
observation in conjunction with the expression for the /-kernels in eq. (80) taken at face value suggests that these 
kernels might not have a well-defined finite Dirichlet limit. However, we evaluated the second order contribution in 
eq. (50). In fact it can be verified that all 2-point kernels have a well-defined finite Dirichlet limit. Further, jl'jj'i' 

also has a well-defined finite Dirichlet limit. The rest of the nine /-kernels on the right hand side of eq. (86) do not 
have a finite Dirichlet limit. But, the sums of the /-kernels listed in eq. (86) have finite Dirichlet limits. The higher 
powers in Ai i2 in the numerator of each of these sums cancel identically to give a well-defined limit as Ai j2 — > oo. This 
seems to be a generic phenomena. Here we list the /-kernels in eq. (86) evaluated in the Dirichlet limit: 



Jd' 2)A («1,«2,«3,«4) 

r(2,2)B/ \ 
P D (Kl,K2,«3,«4) 

r D / (/Cl,K2,K3,«4) 

P D / (Kl,K2,K3,«4) 



«1 



K-2 



«3 



K 4 



sinh K\a sinh n 2 a sinh n 3 a sinh K 4 a 

K\ K 2 K 3 K 4 

sinh K\a tanh n 2 a sinh K 3 a tanh n 4 a ' 



K 4 



1 

sinh Kia sinh n 4 a 4 

Kl K 2 1 

sinh«ia sinh K 2 a 4 



k 2 



K3 



tanh n 2 a tanh n 3 a 
tanh n 3 a tanh K40 



K9 KjO 



+ 



(87a) 
(87b) 

(87c) 
(87d) 



2. Weak coupling limit 



We can similarly evaluate the /-kernels in the weak limit by keeping terms to order A1A2 in the general expressions 
for the /-kernels. But, we find it instructive, and simpler, to start over from eq. (76). We observe that only the 
2-point kernels contribute in the weak limit. We further notice that only the contributions from the free Green's 
function (which is obtained by switching off the couplings Ai, 2 in eq. (43)) contribute in the evaluation of the weak 
limit. These observations allows us to write the 2-point kernels in the weak limit as 



r (mi,m 2 )/ \ ^1 ^2 

(1 2) W l K l' K 2j — 7 7 

K ' h mi! m 2 \ 



Qmi gm 2 e ~Ki\z 2 -zi\ e ~K 2 \zi-z 2 \ 

dz™ 1 dz™ 2 2^ 2^ 



z\=a\ ,z 2 =a 2 



= ^4^r^4^(«i + « 2 ) M e-^+-), (88) 
mi\m 2 \ 2ki2k 2 

where M = mi + m 2 . Note that the derivatives in the above expressions are well defined because they are evaluated 
at a point where z 2 > z\. In particular using these in eq. (86) with the observation that only the 2-point kernels 

contribute in the weak limit we get I^-' 2 ^ A (ki, k 2 , k 3 , k 4 ) = and 

I$ 2)B ( Kll K 2lK3lKi ) = /gJVOci.Ka) = ^ («i + K 3 )*e-«"+">\ (89a) 

/^(M, «a, « 3) «4) = /{JiVfci, K4) = («i + «4) 4 e~ a ^ +K *\ (89b) 

I# 3) (ki, «2, «3, «4) = 4Hw(^^) = YT3T («i + e~ a ^ + «'\ (89c) 

where notice how they differ in the particular dependence on the variables. We again (see comments after eq. (54)) 
point out that the averaging prescription for the evaluation of the /-kernels was not necessary to arrive at the above 
expressions. 
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Using the change of variables introduced in evaluating eq. (55) we can similarly evaluate the L-kernels using eq. (79) 
in the weak limit as 



L {1,2),W ^1,^) mi , m2 , ,g a 



d 



M-l 



I e -«(|fei| + |fc 2 |) 



(90) 



which reproduces the result in eq. (55) for M = 2. 



B. Sinusoidal corrugations: Next-to-leading order 



We shall now specialize to the particular case of sinusoidal corrugations described by eq. (56). Using the Fourier 
transforms in eq. (57) all but one of the k integrals in eq. (85) can be immediately performed. The symmetries listed 
in eq. (83) lead to simplifications in the expressions. As in section III the expressions boil down to two integrals in 
variables k and R. In particular eq. (85) takes the form 



E 



(2,2) 



L X Ly 



■ cos(2k y ) 



h\h\ '' x 
64tt 2 



E 



(3,1) 



L X Ly 



E 



(1,3) 
12 



L X Ly 



= -cos(k y Q ) 



= -cos(fc yo) 



h\h 2 
64tt 2 

hih\ 
64tt 2 



dk I ndn 

oo Jo 

dk RdR 

o JO 



i lW A {n, k+, K, K+) + I^ 2)B (>i-,K, K+, k) 



(91a) 



/oo poo 
dk R 
-oo Jo 



dn 



I^(k, k + ,k,k + ) + / (3 ' 1) (k,k + ,k, k_) + / (3 ' 1) (k,k_,k,k + ) , (91b) 
7 (1 ' 3) (k, k+,k, k+) + / (1,3) (k, k+,k, k-) +/ (1 ' 3) (k, k-,k, k+) , (91c) 



where the /-kernels are given from eq. (86). We use the notations k± and k± introduced after eq. (58). We note the 
factor of 2 in the argument of cosine function of the first term above is a mark of the fourth order in perturbation 
theory. It should be mentioned here that in the above expressions we have omitted finite terms not having any 
dependence in the lateral shift variable y . These terms do not contribute to the lateral force. 



Dirichlet limit 



We presented the expressions for the /-kernels in the Dirichlet limit in eq. (87). Using these in eq. (91) and then 
adding the contributions from the three terms in eq. (74) we can write the total contribution to interaction energy 
due to the presence of corrugations as 



L X Ly 



h x h 2 15 



cos(fc y ) 



240 a 3 a a 4 
where we have introduced the functions 

i poc pC 
-j^ poo poo 

(to) = — sds 

* Jo J-oc 



(hod) 



h 2 ,(1,3) 
— 9 D 



(k a) 



cos 



(2k y ) 1 -^A^\k a) 
Add 



dt- 



.s - _ 



4(2,2) 



dt 



sinh s sinh s + 



sinh 2 s sinh 2 s_ 



tanh s tanh s 
. s 2 

+ 2 



tanh s tanh s+ 



2 2 

- S - S_ 



tanh s sinh.s + sinh.s_ 



(92) 

(93a) 
(93b) 



These functions are plotted in figure 4. As in section III the above functions are suitably normalized such that they 
evaluate to unity when evaluated at to = 0. Verification of this statement involves the following integrals: 



4 3,1) (<>) 



,(2,2) 



(0) 



7r* 
2 



' ds- 



s 2 ds 



sinh s 

A 



tanh s 
.,2 

+ 2 



tanh 2 s sinh 2 s 



( \ 71-2 \ 71-2 

/ 1 2 7T 2 \ f 

U + ^3-) + ( 



2^ 2 - 
63 , 



1. 



sinh s 

Therefore, the next-to-leading-order contribution to the lateral Casimir force in the Dirichlet limit reads 

h\ , h l\ ,(3,i)/, N „ „ x hi h 



Lat.D 



2k a a sin(fc y ) 



5.(0) 
Cas,D 



hi hi 15 
a a A 



+ ^|) Ag' 1) (fcoa)-2cos(fc 2/o) 
cl j d d 



2 /l(2,2)/, n 

A D '(koa) 



(94a) 
(94b) 

(95) 



15 



.(3,1) 




FIG. 4: Dirichlet limit: Plot of (k a) and A^'(k a) versus k a. The dashed curve represents A D ' (koa) which is plotted 

here for reference. 



l(2,2). 



I (1,1)/ 



1.0 




1.0 




FIG. 5: Weak coupling limit: Plot of A { §(k a) and A^(2fe a) versus koa. The dashed curve represents A\J (koa) which is 
plotted here for reference. 

2. Weak coupling limit 

We have already evaluated the L- kernel for the weak coupling case in eq (90), so we can immediately evaluate the 
interaction energy. Alternatively, we can use the expressions for the /-kernels in the weak coupling limit (see eq. (89)) 
into eq. (91) to evaluate the integrals exactly in the spirit of section III and the answer can be expressed as derivatives. 
Using these in (74) we can write the total contribution to the interaction energy in the presence of corrugations in 
the weak limit to be 



E 



(4) _ 
12, W ~ 



E, 



(0) 

Cas : W 



hi /12 3 
a a 2 



cost /„,/„! (4 + 4) A wM -cos(2fc yo)i- -4^(2^/) 



(96) 



where in contrast to the result in the Dirichlet limit we observe the appearance of the same function, though with 
different arguments, as coefficients. This function has been suitably defined as 



4! 84 



1 



-to 



= e ~ t0 E = 

^ ml 

rn—O 



e 4 (t ) 



(97) 



so that it equals unity at to = 0. These are plotted in figure 5. Using eq. (96) in eq. (29) we find that the fourth order 
contribution to the lateral Casimir force in the weak coupling limit equals 



F, 



(4) 



Lat.W 



k a sin(fc 2/o) 



F, 



(0) 



Cas,W 



hi h 2 3 
a a 2 



h 2 h 2 

1 +-|) A^ ) (fc a)-2cos(fc 2/o) 



hi h 2 



A$(2k a) 



(98) 



An exact functional form for the lateral force in the weak coupling limit in the second and fourth order in the per- 
turbative expansion, see eq. (69) and eq. (98), allows us to investigate the approximation involved in the perturbative 
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analysis more closely. We observe that in each perturbative order the leading powers of k^a in the A-function cancels 
the a dependence in the prefactors associated with hi. Thus we conclude that the perturbation analysis is valid for 
koh <C 1. The formal expansions performed in the intermediate steps further restricts our analysis to be valid for 
hi + h 2 < a for arbitrary offset. 



V. PROXIMITY FORCE APPROXIMATION 



As a preliminary check for the results calculated for the lateral force we shall verify that in the proximity force 
limit (koa — ► 0, keeping hi/a fixed,) our results in eqs. (63), (69), (95), and (98), reproduce the standard PFA result 
approximated to appropriate orders in hi /a. For this purpose let us evaluate the lateral Casimir force in the proximity 
force approximation for the sinusoidal corrugations described by eq. (56). We begin by writing the distance between 
the corrugated plates as 

a(y) = a + h 2 sin[fc y] - hi sin[fc (j/ + yo)]- (99) 

For sufficiently small a (in comparison to d = 2ir k^ 1 ) we can treat the plates to be built out of small sections of 
length dy for which the energy is approximately that of the parallel plate geometry. 



1. Dirichlet limit 



Using the expression for the Casimir force between parallel plates in the Dirichlet limit we can thus write 

- 2 1 



dE^ A (a(y)) = -L x dy 



1440 a(y) 3 ' 



(100) 



where we interpret L x dy to be the area of the small section under consideration. Thus the total Casimir energy in 
this approximation, after interpreting L y = hin^^oo Nd, will be 



pPFA 

L X Ly 



1 



1440 a 3 2tt J_ n ^ + h± sin e _ hx sin (0 + ko y Q ) 



Ml d0- 



1440 a 3 2ir 



[1 



(101) 



COS* 



where we have used the substitutions: rsina = h 2 — hi cos koy , rcosa = his'mkoyo, and used the periodicity 
property of the function to eliminate a. We note that r 2 = h\ + h\ — 2h x h 2 cos(k a y a ). 

As was emphasized in the calculations in the earlier sections we have been subtracting the energies that do not 
contribute to the torque (or the lateral force) by treating it as a background. This has not been achieved in the 
above expression in eq. (101) and thus should not be naively compared with our earlier results. But, evaluation of the 
lateral force will let us compare the expressions. Thus, we calculate the lateral Casimir force in the proximity force 
approximation, using eq. (101) in eq. (29), to be 



F EIt A D = 2fc a sin(fcoyo) 



F, 



(0) 



Cas,D 



hi h 2 1 



1 



The basic integral involved in evaluating the above expression is 



for \hi \ + \h 2 \ < a. 



i r 1 1 

2tt J_ v (a -r cosO) ~ Va 2 - r 2 : 



for r < a, 



(102) 



(103) 



which, after the substitution z = cxp(z#), leads to a contour integral on the unit circle in the complex plane. Thus, 
using Cauchy's theorem only those roots of the denominator that are within the unit circle contribute to the integral. 

The PFA result in the leading order for the perturbation parameter is obtained by setting r = in eq. (102). This 
matches exactly with eq. (63) to the leading order in k^a. This verifies the proximity force limit of the leading order 
perturbative result. The proximity force limit of the next-to-lcading-order result in eq. (95) is similarly verified by 
matching it with eq. (102) after evaluating eq. (102) to the fourth order in the perturbation parameter. 
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2. Weak coupling limit 

In the weak coupling limit and in the proximity approximation the equation corresponding to eq. (100) is 



dE^ A (a(y)) = -L x dy 



AiA 2 



32 7T 2 a(y) 

and leads to, after extracting L y , and proceeding in the same manner as in eq. (101), 

AiA 2 1 1 f n ™ 1 AiA 2 1 1 



pPFA 

L X Ly 



de 



32 7T 2 a 2tt J_ n (l + £cos(9) 32 tt 2 a (i - J^)h 



for I /z. x I + \h 2 \ < a. 



(104) 



(105) 



We calculate the lateral Casimir force in the weak coupling limit in the proximity force approximation, using eq. (29), 
to be 



F EIt A w = k oa sin(fc yo) 



P (o) 

Cas.W 



hi h 2 1 

~~ (i-^)i 



for + |/i 2 | < a. 



(106) 



The proximity force limit of the leading- and next-to- leading-order results in the weak coupling limit, eq. (69) and 
eq. (98), can be verified by keeping the appropriate terms in eq. (106). 

VI. ANALYSIS AND DISCUSSION 

In this section we shall exclusively deal either in the Dirichlet limit or in the weak coupling limit. We shall further 
confine ourselves to the case h\ = hi = h. Also, we shall find it convenient to define the following dimensionless 
quantities 



T D ( k a, -, k y 



T w ( k a, -, k y 



5.(0) 

Cas.D 



^ k a sm(k y ) 



F La t,w (koa, %,koyo) 



Cas,W 



k a sin(fc y ) 



(107a) 



(107b) 



We shall use superscripts [m, n] on these dimensionless quantities to denote the order in the power series expansion 
to which the lateral force (not T) have been calculated. In particular m will signify the order in the parameter koa, 
and n will denote the order in h/a. The result in the proximity force approximation should be valid to the first order 
in k a. Thus jfI 1 ' 00 ! will represent the PFA result in our notation. The perturbative results are complete in the koa 
dependence and thus will be denoted by J 7 ! 00 '"! . We use square brackets to avoid the confusion it might lead to with 
the already introduced notation to denote the powers in h^s in the superscripts of the A- functions. In these notations 
our results can be summarized as, (see eqs. (63), (69), (95), (98), (102), and (106),) 



00,4] 



A^ X \koa) + ^ A^'ikoa) - cos(k y ) A [ ^> (k a) 



15 h 2 



,(3,1) 



(2,2), 



•> 1 



00,4] 
w 



e 2 (k a) 



okoa 



2 a 2 
h 2 re 4 (fc a) 

*J n 



[1,00] 



D 



ckoa 



cos(k y a ) 



ei{2koa) 



•> 1 



[l,oo] 



1 



w 



(l-^)t 



(108a) 
(108b) 



where we recall the definition of r after eq. (101). As we mentioned earlier these expressions match to the order [1,4] 
both in the Dirichlet and weak limit. Recall, as justified after eq. (98), that koh is the perturbative parameter. Thus, 
we shall find it convenient to keep koh fixed while discussing the perturbative results. This should be contrasted with 
the case of proximity force approximation where it is more convenient to keep the parameter h/a fixed. All of the 
results in eq. (108) are valid for 2h < a. 



Dirichlet limit 



We plot the lateral force for the Dirichlet limit, J 7 ]^ ' 4 ', with the fractional corrections due to the next-to- leading 
contribution, in figure 6. The plots diverge as koa — > 2koh because the perturbation series breaks down beyond that 
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FIG. 6: Dirichlet limit: The plots on the left show J-^' versus koa for various values of koyo- The dashed curve represents 
Jp°' 2 ' which is plotted here for reference. In each plot the higher values of koh deviate more from the reference. The plots on 
the right show the fractional correction in the next-to-leading order contribution. In each plot the higher values of koh have 
larger corrections. 
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point. The plots for fcoj/o = do not diverge, see eq. (108). We observe that, while keeping fcoyo fixed, the fractional 
corrections due to the next-to- leading order are small for k h -C 1. 

In the absence of an exact answer, for the Dirichlct case, to compare with it is not possible to extract the precise 
error in the perturbative results. In the next section we shall evaluate the corresponding result in the weak coupling 
limit and evaluate the lateral Casimir force non-perturbatively. We show that the error in the lateral force, for the 
weak coupling limit, after the next-to- leading order has been included is small for koh <C 1. Presuming that the same 
would hold for the Dirichlet case, we can expect the error to be sufficiently small for comparison with experimental 
data. (Comparison with experiments requires the corresponding evaluation for the electromagnetic case, which will 
be taken up in a later publication.) 



2. Weak coupling limit - Non-perturbative results 



The plots for the lateral force in the weak coupling limit, JF|^' 4 ', which have not been displayed here, have the 
same qualitative content as the Dirichlct case shown in figure 6. To estimate the error in the perturbative results we 
shall now attempt a non-perturbative evaluation of the lateral Casimir force in the weak coupling limit. 

We begin by pointing out the remarkable observation in [17] that exact results for the Casimir force, in the weak 
coupling limit, can be achieved for specific geometries. This is being extended for a class of geometries, for the scalar 
case, in [19], and another paper will cover dielectrics in electromagnetism, in [20]. Starting from the central formula 
in [19], (alternatively using eq. (38), for the weak coupling limit,) we will have 



Ew 



\\\ 2 

327T 3 



/CO I'OO 
dy / dy' 
-co J — CO 



1 



(y -y') 2 + a(y,y') 2 



where a(y, y') = a + h 2 sm{k y'] — hi sin[fc (y + yo)]- Changing variables, y — y' 
substitutions very similar to those made after eq. (101), we have 



(109) 

y, y + y' — > 2#/fc , and making 



w 



L X Ly 



AxA 2 1 
~32tt 2 2tt 2 



dy 



2tt 



d6 



y 2 + [a — r(k y) cos 9] 2 



(110) 



where r(koy) 2 = h 2 + h 2 — 2hih 2 cos(fcoj/ + &o2/o) 7 which is different from the expression for r in the PFA case, presented 
after eq. (101), in the dependence in the argument of the cosine function. Nevertheless, the 6 integral can still be 
performed using eq. (103) to yield 



L X Ly 



AiA 2 1 
' 32tt 2 tt 



^Re 



(y 



r(k y) 2 



for \hi\ + \h 2 \ < a. 



(Ill) 



Using the above expression in eq. (29) we calculate the lateral force per unit area in the weak coupling limit to be 



-pLat.w A1A2 hi h 2 . .4 1 

T T = ~ 00 2 2 ( fc ° a ) ~ 

L x L y 32ir z a z a a tt 



dt 



Re 



sin(i + k yo) 



[(t + ik Q a) 2 + {k r(t)} 2 p 



for \hi\ + \h 2 \ < a. 



(112) 



Expanding the denominator in the above expression as a power series in kor/(t + ik n a), lets us evaluate the integral 
using the residue theorem. We can thus write, after letting hi = h 2 = h, which leads to the simplification, r(k yo) — 
2/isin[(fc y + fc 2/o)/2], 
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L X Ly 
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32tt 2 a 2 
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-ikoVo 1 

(-i) ^ MX 2 2^+ 2 Y"da 
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(h d Y n+2 1 
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e -koa sin 2n 



/k y .k a\ 



(113) 



which reproduces the leading and next-to-leading results in eqs. (69) and (98) for the case n = and n = 1 respectively. 
This structure justifies our presumption that the perturbative analysis holds as an asymptotic expansion as kgh — > 0. 
Expressing the above sum in a closed form, or evaluation of the integral in eq. (112), would lead to an exact expression 
for the lateral force in the weak coupling limit. We have not been able to achieve either, however, so we shall rely on 
numerical evaluation. 

In terms of the notation introduced in eq. (107) we can write eq. (112) as 



^[00,00] = (koa) 1 
w sin(koyo) tt 



dt 



Re 



sm(t + k y ) 



[{t + ik Q a) 2 + {k r(t)} 2 



for 2h < a. 



(114) 
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FIG. 7: Weak coupling limit: The plots on the left show J r ^°'°° versus koa for various values of fcoi/o at koh = 0.3. The dashed 
curves represents J-~^°' 2 ^ and .T-j^ ' 4 ' which is plotted here for reference. In each plot the exact value is greater than the value 
estimated by J-"j^°' 4 '. The plots on the right show the fractional error in the perturbative result. In each plot the higher values 
of koh have larger errors. 
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FIG. 8: Weak coupling limit: Fractional error in PFA is plotted as the dashed curve. The corresponding error in the perturbative 
result is plotted as the bold curve and is described by the axis on the right. 



The above non-perturbative expression for the lateral Casimir force on corrugated plates, in the weak coupling limit, 
has been plotted in figure 7. For comparison purposes the perturbative results are plotted as dashed curves in figure 7. 
Actually, we found the numerical evaluation of eq. (110), after evaluating the derivative with respect to yo, easier 
than that of eq. (114). We observe that the perturbative results, when the next-to-leading order is included, compares 
with the exact result remarkably well for k^h <C 1 and 2h < a. The fractional error in the lateral force when the 
next-to-leading order contribution is included is displayed on the right column in figure 7. 



3. Validity of PFA 

We shall now compare the PFA result in the weak limit, derived in eq. (106), with the non-perturbative result in 
eq. (112). In figure 8 we plot the fractional error in the PFA result, for fixed koh, and compare it to the fractional 
error in the perturbative results. We note that PFA is a good approximation for koa <C 1 and perturbative results 
are valid for k^h <§C 1. Both the PFA and the perturbative analysis are restricted to 2h < a. The 2h < a restriction 
is necessary for the force being finite, due to non-overlapping of the potentials for arbitrary offset. For special offsets, 
the restriction is weaker, and there is no restriction at all, for zero offset. The precise dependence of the restriction 
on the offset is gotten by referring to eq. (106) and observing that the validity regime is r = 2h sin(/coyo/2) < a. 

A recent debate, see [11-13], involved the comparison of the PFA result and the perturbative result. We earlier 
noted that the PFA limit is obtained by taking the limit koa — ^ while keeping the ratio h/a fixed. In figure 9 we 
plot the lateral force in the PFA limit and compare it with the non-perturbative result for various values of koa. This 
justifies our presumption that the PFA is a good approximation for koa <C 1. We note that the error in the PFA is 
less than 1% for koa <C 1 for arbitrary h/a. For koyo = 7r /4, we observe that the PFA is a very good approximation 
for h ~ a which satisfies 2/isin(fcoJ/o/2) < a. After viewing the plots for various values of koya we note that in general 
the PFA is a good approximation for h ~ a and further beyond for offsets koyo < tt/4. It is, in fact, plausible that 
the PFA holds for large amplitude corrugations for small offsets because the corrugations fit together like fingers in a 
glove. 
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FIG. 9: Weak coupling limit: On the left the PFA is shown as the dashed curve. The exact results for higher values of koa are 
seen to deviate from the dashed curve. On the right the fractional error in PFA is shown to increase with higher values of koa. 



APPENDIX A: DERIVATIVES OF THE GREEN'S FUNCTION 



Evaluation of eq. (47) involves taking derivatives of the Green's function which as we shall see involves evaluating 
a function at a point where it has jump discontinuities in the derivatives. We start by noting that the differential 
equation in eq. (43) can be solved in terms of exponential functions in the different regions described in figure 10. 
Explicitly the solutions are given by the following piecewise-defined functions, with subscripts denoting the regions 
they are defined in, 
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Using the reciprocal symmetry of the Green's function we further have (z, z'\ k) = g^ (z' , z; k), g^'^z, z'\ k) = 
gf\z', z; k), and <?g (2, z'; k) — g?(z', z; k). The above piecewise solution gives the complete Green's function. 

It is of relevance to observe that the Green's function represented above is continuous everywhere in the z, z' domain 
while its first derivatives have simple (jump) discontinuities along the lines, z — z', z — ai.2, and z = ai,2- Since the 
evaluation of /-kernels in eq. (47), and in eq. (80), involves taking derivatives of the Green's function at the lines of 
discontinuity, a definite prescription for the evaluation of the derivatives is called for. For this purpose, we interpret 
the value of a function at the point where it has a jump discontinuity to be the average of all the possible limiting 
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FIG. 10: Description of the regions in the z — z domain on which the piecewise functions contributing to the Green's function 
are defined. 



values. For example, the derivative of the Green's function at the point (z = a%, z' = Oi) in figure 10 takes on four 
different values when approached from the four regions 2, 6, 8, and 9. The value of the derivative is defined as the 
average of these four values. 

We provide the following supporting arguments for the averaging prescription. Firstly, we note that eq. (47) was 
derived from eq (41) after integrating over the complete z — z' domain. The potentials in eq (41) involved derivatives 
of delta functions, and after integration by parts the integrals got contributions from a single point in the z — z' 
domain. Since an integral is a sum we expect it to support the averaging prescription. The second justification we 
provide is in the way of verification for known examples. We made checks by using this prescription to analyze the 
Green's function for the case of a single plate. We do not provide the details of the exercise here. 

Evaluation of the derivatives using the averaging prescription is made convenient when we define the Green's 
function around a point as the average of the suitable piecewise functions. Thus, we evaluate 
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The first derivatives evaluate to be 
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where 3 Z , Z ' is derivative with respect to z, z'. The second derivatives involving two distinct variables evaluate to be 



A 2k 2k 

1 K A l _-na 



\ Z — 0,2, Z — C12 



d z d z ,gW(z,z';K)\ 
d z d z ,gW(z,z';K)\ 



1 n z 



z—ai ,z' —a2 



z—a\,z'—a\ 



z=a 2 ,z' =a 2 



A 2k 
1_k2 
A 2k 

1_k2 
A 2k 



Ai 



2 k 



A 2 



2k 



Ai 



2 k 



Ai 



2 k 



2k 



2 k 



2 k 



Ai _ 2 



2k 



(A5a) 
(A5b) 

(A5c) 



24 



The second derivatives involving the same variables evaluate to 

d 2 J°\z,z'; K ) = K 2 gW(z,z';K) (A6) 

which when compared with eq. (43) tells us that the averaging prescription practically throws away the delta function 
contributions in eq. (43). A more detailed study of this issue is being sought. Since the second derivative returns the 
Green's function back, all higher order derivatives are obtained in terms of the above expressions. 



APPENDIX B: COMPARISON OF OUR LEADING ORDER RESULT WITH REFERENCE [6] 



Here we wish to explicitly compare our results — in the leading order — with those in [6] . In particular we show how 
the Dirichlet limit of our expression for the interaction energy in the leading order [eq. (44) with eq. (51) inserted] 
matches with the results in [6] [eq. (11) with eqs. (13) and (17) there]. To this end we begin from eq. (41) and write 
the interaction energy in the form 
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where we have used the notation introduced in [6]: y = (t,x,y), |y — y'| 2 = (t — t') 2 + {x — x') 2 + {y — y') 2 , and switched 
to Euclidean time by replacing t, t', t — > —it, —it', —it. The Q-kernel introduced above, following [6], is given as 
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where x = (x, y, z) and (x, t, x', t') is given in terms of eq. (42) after switching to Euclidean time as 
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Dirichlet Limit 



In the Dirichlet limit the above expression should correspond to the Q-kernel introduced in [6] [eq. (13) there]. 
Note that in [6] the interaction energy is not isolated from the outset unlike in our eq. (37). This is achieved in [6] 
by identifying and subtracting the contributions to the energy from the single plates. We identify the presence of our 
/-kernel [see eq. (47)] in eq. (B2) which lets us express the Dirichlet limit of the Q-kernel, using eq. (50), in the form 
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This exact form for the Q-kernel in the Dirichlet limit reproduces the result in [6]. 

This exact form for the Q-kernel in the Dirichlet limit suggests that we can consider a class of corrugations for which 
exact results, in the perturbative approximation, might be achievable. Like the exact results achieved in [19, 20] for the 
weak coupling limit, we should be able to explore geometries in the Dirichlet limit in the perturbative approximation 
starting from eq. (Bl). Further, it should also be possible to extend these explorations to the next-to-leading orders. 
And, for very special geometries it might just be possible to transcend the perturbative approximation and get exact 
results in the Dirichlet limit alone. We hope to be able to address some of these explorations in our forthcoming 
publications. 
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